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Current  and  recent  work  in  our  group  on  Josephson  junction  arrays  has 
been  focused  on  two  particular  types  of  dynamical  states,  namely  chaotic 
states  and  “chain”  states.  In  the  following  we  summarize  briefly  the  issues 
addressed  and  the  significant  results  obtained  for  each  type  of  state. 

CHAIN  STATES 

The  dynamical  states  to  which  we  have  given  the  name  chain  states  are 
at  once  both  simple  and  subtle.  The  simplicity,  plus  the  fact  that  this 
type  of  state  is  the  stable  one  in  capacitive  arrays  for  a  broad  range  of 
junction  parameters  and  input  currents,  suggests  that  they  could  be  play  an 
important  role  in  devices  based  on  Josephson-junction  arrays.  The  subtlety 
makes  these  states  interesting  as  examples  of  nonlinear  phenomena  in  their 
own  right. 

The  simplicity  of  the  chain  states  lies  in  their  general  properties  and 
makes  them  easy  to  characterize.  They  have  the  property  that,  given  uni¬ 
form  input  currents  along  two  opposite  edges  of  a  square  array,  array  as  in 
Fig.  1  of  the  accompanying  reprint  of  our  publication  [lj,  the  current  passes 
directly  through  the  junctions  along  the  line  of  the  input  currents,  the  z  di¬ 
rection  in  the  figure.  There  are  no  currents  at  all  in  the  junctions  which  lie 
in  the  transverse,  or  y,  direction.  Further,  the  superconducting  phase  differ¬ 
ence  across  any  longitudinal  junction  t,  meaning  one  which  lies  along 
the  line  of  the  input  current,  has  a  time  dependence  which  is  independent 
of  the  junction,  aside  from  a  translation  in  time.  That  is,  4>i(t)  =  <po(t  ~  ”>»). 


1 


Further,  rt-  is  the  same  for  junctions  separated  from  one  another  in  the 
transverse  directions.  And  finally,  the  same  function  <f>o(t-r)  describes  the 
phase  difference  in  a  single  junction  if  it  has  the  same  parameters  and  is 
driven  with  the  same  input  currents  as  in  the  array.  The  function  ^o(t)  is 
quasiperiodic  in  time  as  opposed  to  being  periodic  or  chaotic.  These  obser¬ 
vations  follow  from  study  of  the  Poincare  plots  of  the  various  junctions  in 
question  or  of  their  power  spectra. 

It  is  easy  to  show  from  analytic  study  of  the  equations  of  motion  for 
the  array  using  the  RSJ  model,  cf.  Eq.  (13)  of  the  print,  that  a  state  of  the 
sort  just  described  is  a  solution  if  Tj  =  njTo  where  n  is  an  integer  and  r0  is 
the  period  of  the  rf  input  current.  What  is  surprising,  and  subtle,  is  the 
stability  of  these  dynamical  states.  For  example,  in  numerical  solutions  of 
the  equations  of  motion,  we  find  that  if  the  array  is  started  from  a  uniform 
state,  meaning  a  state  with  the  same  function  4>o(t)  in  all  longitudinal  junc¬ 
tions,  then  in  the  absence  of  noise  this  state  will  persist  for  an  arbitrarily 
long  time.  If  the  array  is  started  with  random  initial  conditions,  then  it  will 
settle  into  a  chain  state  with  an  apparently  random  but  constant  distribu¬ 
tion  of  r’s,  i.e.,  of  n’s.  In  effect,  the  basin  of  attraction  of  any  given  chain 
state  is  rather  small,  so  that  it  is  extremely  unlikely  the  array  will  wind  up 
in  a  given  chain  state  when  started  from  some  arbitrary  initial  conditions; 
at  the  same  time,  the  array  apparently  always  winds  up  in  some  chain  state. 

Of  course,  in  any  real  array,  there  are  sources  of  noise,  such  as  finite 
temperature  effects,  and  it  is  of  some  importance  to  address  the  question 
of  stability  of  the  chain  states  under  these  circumstances.  From  our  work, 
it  is  clear  that  a  chain  state  with  a  fixed  set  of  n’s  is  not  stable  in  that 
a  burst  of  noise  will  change  the  relative  time  shift  of  the  various  junctions 
in  the  longitudinal  chain.  However,  the  array  as  a  whole  will  return  to  a 
chain  state  with  a  different  set  of  n’s  if  the  noise  is  removed.  If  the  noise 
is  not  removed,  the  system  will  remain  close  to  a  chain  state  in  the  sense 
that  the  currents  in  the  transverse  junctions  will  typically  be  about  the  size 
of  the  noise  currents  while  the  currents  in  the  longitudinal  junctions  v.  ill 
retain  much  of  the  character  they  had  in  the  absence  of  the  noise  as  can  be 
ascertained  by  study  of  the  individual  junctions’  power  spectra.  However, 
the  relative  time  shifts  of  the  phases  of  the  junctions  in  a  chain  will  change 
continually  in  a  random  fashion. 

The  global  properties  of  the  arrays  ,such  as  I-V  characteristics  also  merit 
some  discussion,  especially  when  there  is  noise  present.  The  noise  and  the 
fact  that  many  junctions  contribute  to  the  global  properties  have  the  effect 
of  causing  the  long-time  state  of  the  system  to  show  very  little  dependence 
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on  the  initial  state.  Moreover,  the  voltage  across  the  array  mimics  closely 
the  input  current  in  that  if  the  input  is  Ijc  +  /,/  sin(u>oO>  then  the  output 
voltage  is  of  the  form  Vjc  +  VTj  sin(u>0t  +  a),  i.e.,  it  is  ohmic  aside  from  a 
time  shift.  That  is  rather  remarkable  because  the  response  of  an  individual 
junction  is  definitely  not  linear  as  may  be  seen  from  the  power  spectra  of 
the  junctions  in  the  absence  of  noise  in  the  chain  state,  cf.  fig.  4(b)  of 
Ref  1.  Typically,  there  are  two  incommensurate  frequencies  of  response 
which  show  up  clearly  along  with  their  differences,  harmonics,  etc.,  leading 
to  aperiodic  (quasiperiodic)  behavior  of  individual  junctions.  The  effect  of 
having  many  junctions  in  a  ehain  contributing  to  the  array  voltage,  with 
random  relative  time  shifts,  is  to  remove,  in  the  limit  of  large  arrays,  the 
contributions  from  frequencies  other  than  the  driving  frequency.  That  is,  if 
one  looks  at  the  Fourier  decomposition  of  the  voltage  across  the  array,  it 
shows,  aside  from  a  certain  amount  of  background  noise,  a  large  peak  at  the 
driving  frequency  <j,  while  the  other  peaks,  so  evident  in  the  spectrum  of  a 
single  junction,  are  much  reduced,  apparently  varying  as  1/JV,  where  N  is 
the  linear  size  of  the  array  in  the  x  direction.  Consequently,  a  large  noisy 
(finite  temperature)  array  may  be  expected  to  produce  faithfully  essentially 
ohmic  behavior  with  a  large  voltage  signal  (proportional  to  N)  at  the  same 
frequency  as  the  applied  rf  field  even  though  the  response  of  individual 
junctions  is  dramatically  different. 

SPA  TIO  TEMP  ORA  LL  Y  CHAOTIC  STATES 

Over  the  past  few  years  it  has  been  known  that  a  harmonically  driven 
single  capacitive  Josephson  Junction  exhibits  chaotic  behavior.  Within  the 
Resistively  Shunted  Junction  (RSJ)  model  a  square  array  of  these  junctions 
is  known  to  exhibit  variety  of  dynamical  states  including  chaotic  states  when 
harmonically  driven  by  uniform  input  currents. 

We  have  studied  in  a  recent  preprint[2]  aspects  of  spatiotemporal  chaos 
in  Josephson  Junction  arrays  and  discussed  the  dynamical  behavior  of  both 
spatially  averaged  (  collective  )  and  local  quantities  in  the  chaotic  states. 

One  of  us  (C.J)  in  a  collaborative  effort  a  few  years  ago  [3]  has  argued 
heuristically  in  the  context  of  coupled  lattice  models  that  collectively  chaotic 
states,  e.g.,  ones  in  which  some  fourier  component  is  chaotic,  do  not  occur 
generically  in  systems  with  short  ranged  interactions  in  the  presence  of  ex¬ 
ternal  noise  (  e.g.  Langevin  noise  ).  We  have  confirmed  numerically  that 
collectively  chaotic  states  do  not  occur  in  JJA  described  by  the  the  RSJ 
model  in  agreement  with  the  heuristic  argument.  It  should  be  noted  that 
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the  RSJ  model  is  different  from  the  coupled  lattice  map  models  studied  ear¬ 
lier  in  that  current  conservation  allows  instantaneous  response  of  the  entire 
system  to  local  perturbations.  The  route  to  chaos  is  similar  to  that  of  a 
single  junction  (  analysis  of  the  dynamical  behavior  near  transition  regions 
is  complicated  ).  In  the  chaotic  states  spatial  averages  of  the  quantities  such 
as  voltage  difference  and  josephson  currents  along  parallel  or  perpendicular 
direction  of  the  input  current  do  not  exhibit  chaotic  behavior  in  the  large 
system  size  limit.  (See  Fig.  1)  These  variables  are  periodic  in  time  with  a 
period  set  by  the  external  driving  currents  (tr/sin(u>t)).  The  spread  in  the 
voltage  scales  like  1  /(JV(L))1/3,  where  N(L)  =  L(L  —  1)  is  the  number  of 
junctions  along  the  direction  of  the  input  current  in  a  L  X  L  lattice.  Sim¬ 
ilar  behavior  is  seen  for  other  spatially  averaged  quantities.  In  the  chaotic 
states  the  largest  Lyapunov  exponent  A (L)  has  an  asymptotic  value  Ao  for 
large  L.  For  small  L  the  system  behaves  like  a  linear  array  of  junctions  with 
zero  transverse  currents  and  the  value  of  A  is  close  to  the  value  for  a  single 
junction  The  largest  Lyapunov  exponent  is  shown  as  a  function  of  the  size 
of  the  system  in  Fig.l.  The  cross  over  from  small  to  large  system  behavior 
indicates  a  length  which  is  of  the  order  of  the  length  obtained  from  the 
small  r  decay  of  the  correlation  function  G(r)  that  we  have  also  computed. 

Second,  we  have  studied  the  probability  distribution  functions  (PDF)  for 
the  local  voltage  variables  and  find  that  they  have  exponential  tails  in  the 
chaotic  states.  Such  probability  distributions  for  local  variables  are  observed 
in  fluid  systems  [4].  A  model  Fokker-Planck  equation  has  been  proposed  as 
a  plausible  mechanism  for  these  exponential  tails  [5].  We  have  shown  that 
the  observed  exponential  tails  of  PDF  can  be  understood  by  a  similar  model 
that  describes  the  mixing  of  the  local  voltage  variable  in  the  presence  of 
the  supercurrent  fluctuations.  In  Fig.2  we  show  the  probability  distribution 
function  for  the  voltage  fluctuations  of  a  transverse  (y  dir)  junction  in  the 
chaotic  state.  The  data  convincingly  shows  the  presence  of  the  exponential 
tails  in  the  PDF.  In  contrast  the  PDF  for  local  voltage  variable  fluctuations 
in  a  periodic  state  with  external  noise  (  <r  ^  0  )  is  gaussian  (  see  Fig.2b  ). 
This  can  be  shown  by  linearizing  (2)  around  the  periodic  state  and  writing 
the  equation  of  motion  for  a  local  voltage  variable  in  the  form  of  a  Langevin 
equation. 
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1  Figure  captions 

Fig.l.  Largest  Lyapunov  exponent  as  a  function  of  the  size  of  the  system  for 
id e  =  0.69,  s,/  =  1.154,  /3C  =  4.0,  and  u  =  0.1  X  2x.  Notice  that  it  increases 
for  small  L  and  approaches  an  asymptotic  value  for  large  L. 

Fig.2.  a)  Probability  distribution  for  the  fluctuations  in  the  voltage  of  a 
transverse  junction  for  the  same  parameters  as  in  Fig.l.  b)  shows  the  PDF 
for  the  same  variable  in  a  uniform  nonchaotic  state  in  the  presence  of  noise. 
The  parameters  are  i<u  =  0.4, *,/  =  1.0, 0e  =  4.0,  and  u  =  0.1  x  2x.  Noise 
value  (cr)  =  0.1 
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We  have  studied  the  dynamics  of  regular  two-dimensional  Josephson-junction  at  rays  subjected  to  elec¬ 
tromagnetic  radiation  at  frequencies  comparable  to  the  individual  junction’s  characteristic  frequency. 
The  junctions  are  described  using  the  resistively-shunted-junction  model  including  capacitance  with  the 
plasma  frequency  also  comparable  to  the  characteristic  frequency.  The  dynamical  behavior  fails  into 
several  different  general  classes,  namely,  periodic,  quasiperiodic,  and  chaotic,  depending  on  the  particu¬ 
lar  characteristics  of  the  junctions,  the  input  currents,  and  the  amplitude  and  frequency  of  the  radiation. 
Detailed  examples  of  each  of  these  types  of  behavior  are  given.  Current-voltage  characteristics  are  ex¬ 
amined  and  related  to  the  dynamical  behavior  and  the  junction’s  properties.  The  effect  of  finite  tempera¬ 
tures,  included  by  means  of  a  Langevin  noise  current,  is  also  discussed,  as  is  the  stability  of  various  types 
of  dynamical  states. 


I.  INTRODUCTION 

The  dynamical  behavior  of  Josephson-junction  arrays 
has  been  studied  in  considerable  detail,  both  experimen¬ 
tally1-4  and  theoretically.5  - 10  Emphasis  in  recent 
theoretical  work  has  centered  on  the  dynamics  in  the 
presence  of  an  applied  static  magnetic  field  and  in  the 
limit  of  zero-junction  capacitance  so  that  the  junctions 
are  overdamped.  An  exception,  and  the  work  most  close¬ 
ly  related  to  our  research,  is  Ref.  10  in  which  the  types  of 
dynamical  states  and  their  stability  are  studied  for  a 
linear  series  array  of  junctions  with  capacitance. 

In  this  paper  we  present  some  results  from  an  extensive 
and  continuing  study  of  the  dynamics  of  regular  two- 
dimensional  Josephson-junction  arrays  on  a  square  lattice 
in  the  presence  of  an  incident  electromagnetic  (rf)  field 
and  applied  dc  bias  currents.  A  finite  junction  capaci¬ 
tance  is  included.  Our  intention  is  to  give  a  broad  survey 
of  the  possible  behaviors  of  the  arrays;  more  detailed  ex¬ 
ploration  of  particular  types  of  states  will  be  given  sepa¬ 
rately.  Our  study  is  motivated,  in  part,  by  the  possibili¬ 
ty1  11  of  constructing  detectors  or  generators  of  mi¬ 
crowave  radiation  from  arrays.  A  second,  more  theoreti¬ 
cal,  motivation  is  the  possibility  of  observing  dynamical 
e  g.,  chaotic  and  quasiperiodic,  behavior  of  an  unusual 
sort  in  these  nonlinear,  continuous  time,  coupled  systems. 
We  focus  on  the  character  of  the  locally  stable  modes  of 
oscillation  given  different  input  currents  and  incident  rf 
fields  for  arrays  with  various  characteristic  parameters. 
Numerous  similar  investigations  have  been  performed  for 
a  single  capacitive  Josephson  junction.12, 13  Also,  for  use 
as  detectors,  an  important  measure  of  the  response  of  the 
system  is  the  mean  voltage  across  the  array  and  its  varia¬ 
tion  with  changes  in  the  rf  field  incident  radiation.  Con¬ 
sequently,  we  also  present  typical  l-V  curves  for  the  ar¬ 
rays  and  relate  the  characteristics  of  these  curves  to  the 
dynamical  states. 

Section  II  of  the  paper  contains  descriptions  of  the 
model  employed  for  the  array,  the  numerical  methods 
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used  for  its  analysis,  and  the  quantities  studied;  Sec.  Ill  is 
devoted  to  presentation  of  the  numerical  results  while 
Sec.  IV  is  given  to  a  linear  stability  analysis  of  certain 
dynamical  modes;  a  discussion  and  summary  compose 
Sec.  V. 


II.  MODEL  AND  METHODS  OF  ANALYSIS 

We  describe  a  Josephson  junction  using  the  resistively- 
shunted-junction  (RSJ)  model14  including  the  junction  ca¬ 
pacitance  and  noise.  The  essence  of  this  model  is  that  the 
current  through  a  junction  between  two  superconductors 
is  written  as  a  sum  of  several  kinds  of  currents  in  parallel. 
These  include  a  displacement  current  ID  through  a  ca¬ 
pacitance  C;  a  normal  or  quasiparticle  current  /iV 
through  a  resistance  R;  a  supercurrent  (the  Josephson 
current)  /s;  and  a  fluctuation  or  noise  current  IF.  Thus, 
the  total  current  in  the  junction  is 

/=/fl+/;V+/5  +  /F  .  (1) 


The  displacement  and  normal  currents  may  be  ex¬ 
pressed  in  terms  of  the  voltage  V  across  the  junction, 

,dV 


ID=C 


and 


dt 


(2) 


/v  =  — 


The  supercurrent  is  expressed  as  a  function  of  the  super¬ 
conducting  phase  difference  6  across  the  junction 

Is=Icsind>  ,  (3) 


where  fc  is  the  critical  current.  Further,  the  relation  be¬ 
tween  the  voltage  and  phase  is 


2e  dt 


(4) 
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Finite  temperatures  can  be  introduced  by  means  of  1F 
which  is  taken  to  be  a  Langevin  noise  current  represent¬ 
ing  random  quasiparticle  currents  tunneling  through  the 
junction.  This  current  is  given  the  properties 

</f(t)>=  0 


and  (5) 

(IF(t  +t')IF(t))  =  ^L8(f)  , 

R 

where  the  brackets  (  •  •  •  )  denote  an  ensemble  average; 
T  is  the  temperature  and  k  is  Boltzmann’s  constant. 

Combining  these  relations  and  dividing  by  the  critical 
current,  one  finds  that  the  current  in  the  junction  can  be 
written  as 


ca 

d 2<b  , 

fi 

2  elc 

2eRIc 

d<j> 

dt 


+sin$+/f(f) 


(6) 


where  iF  =IF/IC.  Notice  that,  from  Eq.  (5), 

(iF(t  +t')iF(t)}  —aS(t')  (7) 

with  a  =  2kT/RI'. 

Introduce  the  junction’s  characteristic  frequency 
(oc  =  2eRIc /fi  and  plasma  frequency  a>p  =  \/ 2elc /Cfi\ 
also,  introduce  a  dimension-free  time  r  —  opt.  Then  the 
current  may  be  written  as 

i  =  +0c  W2~r~  +  sind  +  /F  ,  (8) 

dr  dr 

where  0C  =icoc  /a>p  )2.  Equation  (8)  is  the  basic  equation 
of  motion  for  the  RSJ  model  of  a  capacitive  junction. 

Given  an  array  of  junctions  meeting  at  a  set  of  super¬ 
conducting  nodes  k ,  one  can  write  the  RSJ  equation  for 
an  individual  junction  in  terms  of  the  difference  of  the 


phases  of  the  superconducting  order  parameters  on  the 
individual  nodes.  We  assume  a  uniform  array,  meaning 
that  all  junctions  have  the  same  parameters  cop,  o>c,  and 
Ic.  The  equation  for  the  current  in  the  junction  between 
neighboring  nodes  k  and  /  is 

•  _  ^  d2<t>i  ,  a- 1/2  d d^' 

*"  dr2  dr2  Pc  dr  dr 

+%m(<j>k—<^i)+(iF)ki  .  (9) 

We  further  assume  the  nodes  are  sufficiently  small  that 
the  phase  of  a  node  is  spatially  uniform  and  that  charging 
effects  on  a  node  are  negligible  so  that  the  sum  of  the 
currents  entering  a  given  node  is  zero.  This  condition 
determines  how  the  equations  of  motion  of  the  various 
junctions  are  coupled, 

2  *ki~~Ueh  >  (10) 

ie.vik\ 

where  the  sum  is  over  those  neighbors  l  of  k  which  are 
coupled  to  it  by  junctions,  and  (ie)k  is  the  input  current, 
in  units  of  Ic,  to  the  kth  node  from  external  sources. 

Equations  (9)  and  (10)  are  the  equations  of  motion  for 
the  array.  Given  some  initial  conditions  on  the  phases 
and  their  time  derivatives,  a  set  of  input  currents,  and  a 
formulation  of  the  fluctuation  currents,  they  present  a 
well-defined  dynamical  problem  involving  the  solution  of 
a  set  of  coupled,  second-order,  ordinary,  nonlinear 
differential  equations.  They  may  be  expressed  in  matrix 
form15  as 

G  =-sin <f>—iF—ie  ,  (11) 

dr  dr 

where  G  is  a  matrix  and  4>,  sin<£,  iF,  and  it  represent,  re¬ 
spectively,  vectors  comprising  the  <f>k' s,  sink’s,  (iF)k’s, 
and  ( ie  )k ’s.  The  elements  of  G  are 


nk  if  l  —k  and  nk  is  the  number  of  junctions  entering  node  k, 

Gk/~  —  1  if  k=£\  and  k  and  /  are  joined  by  a  junction,  (12) 

0  otherwise  . 


For  any  given  G,  the  inverse  can  be  obtained  numerically. 
There  is  however,  one  technical  point  that  should  be 
mentioned.  If  there  are  N  nodes,  then  there  are  N  equa¬ 
tions.  The  matrix  G  for  this  set  of  equations  possesses  no 
inverse,  a  fact  that  can  be  related  to  overall  current  con¬ 
servation,  meaning  that  the  net  current  entering  the  array 
must  he  zero,  and  to  the  fact  that  the  system  is  invariant 
under  a  shift  of  the  phases  of  all  of  the  nodes.  We  can  ar¬ 
bitrarily  set  the  phase  of  one  node  equal  to  a  constant, 
thereby  reducing  the  number  of  independent  variables 
and  the  number  of  equations  by  one.  The  matrix  G  for 
the  remaining  iV  -  1  equations  possesses  an  inverse  which 
is  easy  to  compute.  Given  this  matrix,  we  recast  the 
equations  of  motion  as 


d!±+p-w2d±=G-'(-sin<t>-iF+ie)  .  (13) 

dr  dr 

This  form  is  more  amenable  to  efficient  numerical  solu¬ 
tion  than  the  original  one. 

We  turn  now  to  the  inclusion  of  incident  electromag¬ 
netic  fields.  A  static  applied  magnetic  induction  and  in¬ 
cident  electromagnetic  radiation  may  be  described  by 
means  of  vector  potentials  As  and  Arf.  The  coupling  to 
the  array  is  introduced  by  replacing  the  phase  difference 
in  the  Josephson  current  <f>k—  d/  with  4>k—<j>i—  A  (k,l), 
where 

Alk,l)=l f  f  dl-  A, 

<P0  Jk 


(14) 
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d>0  is  the  flux  quantum  he  /2e.  In  the  case  of  the  radia¬ 
tion  field,  we  are  interested  in  microwave  frequencies. 
The  corresponding  wavelength,  or  distance  over  which 
the  vector  potential  varies  appreciably,  is  on  the  order  of 
a  centimeter.  We  will  suppose  that  this  is  much  larger 
than  the  size  of  the  array  and  so  shall  approximate  the 
vector  potential  at  all  points  on  the  array  by  that  at  its 
center,  in  effect  making  the  dipole  approximation.  Thus, 
we  employ,  for  a  monochromatic  wave  of  frequency  a)r 
such  as  may  be  obtained  in  a  resonant  cavity, 

Arf(r)=  AoSimu,./  =  Aosin0y2T  ,  (15) 

where  [3r  =  ( tor  /iop  )2.  For  this  case, 

.4  (*,/)=  It  \  .(T  -r  )  .  (16) 

<t>o 


<•  (■*> 


(1,»)  (2.M) 


OH) 


n 


(HU) 


i,(3) 

i.(2) 

<1J) 

03) 

(3,3)  ’  *  * 

(K3) 

(1Z) 

02) 

ozT 

dW 

•  •  • 

k _ 1 

(1.1) 

OD 

(3,1) 

(HD 

X 


Given  an  ordered  array  such  that  the  nodes  are  on  the 
sites  of  a  regular  lattice,  and  with  junctions  only  between 
nearest-neighbor  sites,  the  A(k,l)  are  the  same  for  all 
junctions  oriented  in  a  particular  direction. 

Because  of  the  special  character  of  A  (k,l)  in  the  dipole 
approximation,  a  change  of  variables  can  be  made  to  re¬ 
move  all  reference  to  the  vector  potential  from  the  teims 
representing  the  Josephson  currents.  Define 

<b'k  =6k  +  Ao-rksin(0'/2r)  .  (17) 

***() 


In  terms  of  this  variable,  the  current  in  a  junction  be¬ 
tween  nodes  kj  may  be  written,  with  a  shift  in  the  origin 
of  the  time  variable,  as 


'*/  = 


dli'k 

dr2 


+  Pc 


1/2 


dr 


di>\ 

77 


+  sin((^  —<b])  +  iF 

-  A0-<r*  -r, ) Vfi+PrPc  'sinU3'/2r)  . 

*7) 


(18) 


An  applied  uniform  static  magnetic  induction  B,  can 
be  included  by  adding  a  simple  vector  potential,  e.g., 
A,  =(B,Xr)/2.  Then  the  A(k,lYs  will  depend  on  the 
positions  of  the  sites  (k,l)  and  a  transformation  of  the 
phases  will  not  result  in  simplification  of  the  equations  of 
motion. 

In  the  work  described  in  Sec.  Ill,  we  employ  M  X  Af  ar¬ 
rays  on  a  square  lattice,  as  shown  schematically  in  Fig.  1, 
with  a  distance  a  between  nodes.  For  definiteness,  let  the 
array  lie  in  the  z  =0  plane  with  the  x  direction  along  the 
side  of  the  lattice  of  length  Na  and  the  Y  direction  along 
the  side  of  length  \fa.  We  identify  a  node  with  two  in¬ 
dices  (/',/>,  where  i  =  1,2, . . .  ,A'and  j  =  1,2, ...  .Af.  We 
apply  input  currents  along  the  opposite  edges  /  =  1  and 
I  -  V  These  currents  may  be  functions  of  j.  Total 
current  conservation  is  implemented  by  extracting  the 
same  total  current  from  one  edge  as  is  injected  at  the  op¬ 
posite  edge.  Along  the  two  lateral  edges  of  the  array, 
given  by  j  —  1  to  j  =  .Vf,  we  allow  no  current  to  enter  or 
leave  the  system. 

Because  of  the  orientation  of  the  array,  only  the  x  and 
v  components  of  the  vector  potential  will  enter  the  equa- 


FIG.  1.  Schematic  diagram  of  a  NXM  array.  The  nodes  are 
represented  by  discs  with  the  i  and  j  values  of  a  node  shown  in 
parentheses.  The  junctions  are  represented  by  solid  lines.  Input 
currents  are  also  shown. 


tions  of  motion.  For  an  applied  rf  field,  and  in  the  dipole 
approximation,  we  choose 

Arf= /40(cosax+sinay)sin(j3j/2f)  ,  (19) 

a  is  the  angle  between  Arf  in  the  x-y  plane  and  the  x 
direction.  Then,  from  Eq.  (16), 

Ari(i,j\i  Ancosa sin(0'/2r)  (20) 


and 

At((i,j;i,j  + 1  )=  ^0sina sin(0}/2r)  •  (21) 

When  the  transformation  Eq.  (17)  is  made  to  the  vari¬ 
ables  <b\  we  find  that  this  vector  potential  is  equivalent  to 
an  ac  input  current  along  the  edges  of  the  array  given  by 

I* rf <  l.;>= - - - Vp\  +PrP7  ‘sin (P\/2T) 

P0 

=  —  trf(  A,y)  =  i„sin(j8‘/2T)  (22) 

and 

2ir  A  0a  sina  - ;  ... 

,rf(  /,  1 )  = - - - VP\  +PrP7  ‘sin (PV2t) 

=  —iTfU,M)  =  irysm(0y2T)  .  (23) 

These  equations  serve  to  define  in  and  iv.  Note  that 
both  of  these  contributions  apply  at  the  corners  of  the  ar¬ 
ray. 

The  total  input  current  to  the  array  we  take  as  the  rf 
currents  irf  of  Eqs.  (22)  and  (23)  and  dc  currents  ie  intro¬ 
duced  in  Eq.  (10).  Hence,  the  “external”  currents  appear¬ 
ing  in  the  equations  of  motion,  Eq.  (11)  or  Eq.  (13),  be¬ 
come 

i'.=ie+iv{  (24) 

when  <b  is  replaced  by  <b' . 
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III.  CALCULATIONS  AND  RESULTS 


A.  Zero  temperature:  />  =  0 


Our  numerical  procedure  is  to  integrate  Eq.  (13)  using 
a  fourth-order  Runge-Kutta  method.  This  general  tech¬ 
nique,  but  without  including  the  junctions’  capacitance, 
has  been  used  by  numerous  other  workers6-9  to  solve  for 
the  behavior  of  arrays.  The  step  size  in  time  is  typically 
chosen  to  be  around  0.02  of  the  smallest  of  the  three 
periods  2v/o)p,  2n/coc,  and  2 n/cor.  Smaller  step  sizes  are 
frequently  employed  to  check  the  accuracy  of  the  results. 
Runs  lasting  anywhere  from  10  to  several  thousand  such 
periods  are  generally  sufficient  to  determine  the  character 
of  the  asymptotic  behavior  of  the  array  and  to  obtain  reli¬ 
ably  such  properties  as  phase-space  trajectories  and 
Fourier  transforms  of  the  voltage  or  supercurrent  across 
individual  junctions  or  the  entire  array.  The  array  size 
used  in  calculations  is  routinely  8X8,  but  arrays  of  size, 
e.g.,  32X32  or  16X64  have  been  analyzed  to  examine 
specific  features  and  to  check  for  size  dependence.  The 
properties  discussed  in  this  paper  tend  to  be  quite  insensi¬ 
tive  to  the  size  once  it  is  8X8  or  larger.  In  practice,  we 
measure  times  in  units  of  the  period  of  the  applied  rf 
field;  that  is,  we  employ  the  time  variable 

<°r 


Also,  we  express  time  derivatives  using  this  unit  of  time, 
defining 


dd>'  _  2rr  dd>' 
d  r,  ior  dt 


(26) 


Similarly,  frequencies  are  measured  in  units  of  wr. 

There  are  numerous  parameters  entering  the  equations 
of  motion.  These  include  /?,,  A(l,  a,  and  the  various 
dc  input  currents.  We  have  focused  attention  on  the  re¬ 
gime  where  the  three  frequencies  are  comparable  (/?f  —  1 
and  Pr  ~  1 )  and  where  the  dc  input  currents  and  effective 
rf  input  current  are  comparable  to  the  critical  current. 
i,  ~  1  and  irx  ~~  1 ;  A„  is  taken  parallel  to  the  x  direction  so 
that  ify=0.  We  first  studied  the  dynamical  behavior  in 
the  absence  of  noise  currents  (T=  0)  and  then  added 
noise  to  determine  whether  any  significant  changes  in  the 
dynamical  states  occur  over  and  above  simple  fluctua¬ 
tions  superposed  on  the  T=  0  motion.  In  the  following 
we  first  present  representative  results  with  applied  dc 
currents  and  effective  rf  currents  at  zero  temperature 
b  ~0l  with  no  static  applied  field.  We  then  address  the 
effect  of  finite  temperatures. 

We  have  evaluated  and  analyzed  both  detailed  and  glo¬ 
bal  dynamical  properties  of  arrays.  These  include  the 
'ime-dependent  supercurrents  and  normal  currents,  or 
■  nltuges,  in  individual  junctions  and  their  time  averages; 
'ime-dependent  voltages  across  the  entire  array  and  their 
'ime  averages  'yielding  I-V  characteristics!;  and  the 
-"Mtiguration  and  motions  of  vortices.  Also,  in  the  case 
’l  chaotic  states,  the  character  of  the  chaos  and  its  onset 
base  been  studied  in  some  detail.  In  this  paper  we  do  not 
discuss  the  motion  of  the  vortices  or  the  properties  of  the 
- -antic  states. 


We  describe  the  rf  field  in  terms  of  the  equivalent  input 
currents  i„.  As  stated  earlier,  we  suppose  that  iv  —  0, 
meaning  that  the  array  is  oriented  with  the  x  direction 
parallel  to  the  vector  potential,  or  the  electric  field,  of  the 
incident  radiation.  Finally,  for  the  most  part,  we  take  the 
dc  input  current  in  the  nodes  (1  ,j)  to  be  uniform  along 
the  y  direction  (independent  of  j)  and  given  by  the  single 
number  ie.  Then  the  parameters  of  the  system  in  the  ab¬ 
sence  of  noise  are  ie,  in,  pc,  and  Pr.  For  various  sets  of 
these  parameters  we  solve  for  the  motion  starting  from 
random  initial  conditions. 

For  many  parameter  sets,  the  junctions’  motions 
decouple  in  the  sense  that,  at  long  enough  times,  there 
are  no  currents  in  any  junctions  oriented  along  the  y,  or 
transverse  direction.  Thus,  there  are  simply  independent 
lines  of  current  along  the  x,  or  longitudinal  direction. 
For  this  type  of  motion,  one  need  not  study  a  two- 
dimensional  array  to  determine  the  dynamical  properties. 
The  shapes  and  locations  of  the  regions  in  parameter 
space  where  decoupled  states  appear  are  usually  rather  ir¬ 
regular.  For  example,  Fig.  2(a)  shows  points  on  a  grid  in 
ie-irx  space  where  the  motion  did  not  decouple  in  com¬ 
puter  runs  for  times  up  to  2000ir/cor  with  Pc  =  1  and 
P,  =  1.  In  these  runs,  ie  and  in  were  set  equal  to  multi¬ 
ples  of  0.05  from  0.0  and  2.0  and  the  system  was  started 
from  a  random  configuration  of  <t>' s  and  <j  "s.  It  should  be 
emphasized  that,  for  some  of  the  points  at  which  the 
motion  did  not  decouple  in  our  runs,  it  would  have 
decoupled  if  we  had  run  a  (much)  longer  time.  Also,  the 
amount  of  time  necessary  for  the  system  to  settle  into  a 
decoupled  state  depends,  in  some  cases  quite  strongly,  on 
the  initial  conditions.  Indeed,  there  are  instances,  as  in 
the  case  of  a  single  junction,1’  where  the  character  of  the 
long-time  state  depends  on  the  initial  conditions  because 
there  are  numerous  locally  stable  states  of  motion.  This 
behavior  tends  to  be  especially  true,  not  surprisingly,  if 
the  parameters  are  such  that  the  system  is  close  to  a 
boundary  between  regions  of  predominantly  coupled  and 
predominantly  decoupled  states.  Hence,  the  figure 
should  be  interpreted  as  meaning  that,  where  there  is  a 
point,  the  system  is  likely  to  remain  in  a  coupled  state  for 
a  long  time  if  started  with  random  initial  conditions; 
there  may  be,  but  need  not  be,  one  or  more  locally  stable 
decoupled  states  for  the  same  parameters.  Similarly, 
where  there  is  no  point  in  the  figure,  that  means  the  sys¬ 
tem  is  likely  to  remain  in  a  decoupled  state  at  long  times, 
although  there  may  be  locally  stable  coupled  states. 

The  regions  where  the  motion  decouples  depend 
strongly  on  P,  and  Pr.  General  trends  are  somewhat 
elusive.  One  might  expect  that,  if  Pc  -  1  so  that  ioc  -  cop, 
then  the  long-time  dynamical  state  is  more  likely  to  be 
decoupled  if  u>r  is  either  much  larger  or  much  smaller 
than  the  other  frequencies,  i.e.,  Pr  »  1  or  Pr  «  1.  This 
expectation  is  borne  out  to  some  extent.  For  example,  us¬ 
ing  the  same  grid  as  employed  to  obtain  Fig.  2(a)  we 
found  a  small  number — about  25 — of  coupled  states  at 
long  times  with  p(  =4  and  with  either  p,  =  25  or  Pr  = 

In  these  two  cases  the  locations  of  the  coupled  states  are 
distributed  quite  differently.  For  Pr~  25  they  are  widely 
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spread  in  <,,  but  all  close  to  <,=0.7,  whereas  for0r  =  ±, 
they  are  dispersed  at  relatively  large  in  and  small  <,,  im¬ 
plying  that  there  is  no  simple  rule  for  the  probable  loca¬ 
tions  of  coupled  or  decoupled  states. 

Another  trend  we  have  found  is  the  following:  if  the 
capacitance,  or  0,,  is  small,  so  that  the  displacement 


°5|_ 

.  (Cl 

f)  0  —  .  .  .  '■ . -  -  - - — 

no  os  io  is  :.o 

i. 

FIG.  2.  Points  on  a  grid  with  spacing  0.05  in  the  space  of 
and  t„  at  which  coupled  long-time  dynamical  states  are  found 
for  'ai  p  ~P.  -  I,  <bi  P,  -  I  and  pr  ~ 25 ,  and  lei  P,  =4  and 
P,  rr/25.  Arrays  measuring  8X8  were  used  to  obtain  these 
results  as  well  as  those  in  all  of  the  figures  that  follow. 


current  is  not  important  and  the  system  is  overdamped, 
then  too  the  dynamical  state  will  be  decoupled.  Then,  as 
Pc  increases,  more  coupled  states  can  be  expected.  For 
example,  in  Fig.  2(b),  we  show  the  regions  of  coupled 
states  for  fic  =  1  and  f$r  =  25.  In  contrast  to  the  results  for 
the  same  / 3r  and  Pc  =4,  we  find  a  large  region  where  cou¬ 
pled  states  appear  at  long  times.  For  further  contrast,  we 
show  in  Fig.  2(c)  the  behavior  for  pc  =  4  and  /?, =17^/25 
which  demonstrates  that  the  distribution  of  coupled 
states  can  be  quite  complex. 

To  summarize,  all  cases  shown  give  decoupled  states  if 
the  input  currents  /,  and  i„  are  small  compared  to  unity 
and  also  if  <',  is  significantly  larger  than  unity.  Our  gen¬ 
eral  findings  are  that,  when  started  from  some  random 
configuration,  the  arrays  show  long-time  behavior  in 
which  the  currents  either  travel  straight  through  the  ar¬ 
ray,  so-called  decoupled  states,  or  travel  in  irregular 
paths  so  that  there  are  sizable  currents  in  the  transverse 
junctions,  so-called  coupled  states.  If  fic,  Pr,  <,,  and  in 
are  of  order  unity,  it  is  difficult  to  predict,  in  general, 
which  type  of  behavior  will  be  observed  in  the  absence  of 
detailed  numerical  study  of  that  particular  case.  Howev¬ 
er,  in  some  instances  the  behavior  is  already  determined 
from  knowledge  of  the  dynamics  of  a  single  junction  hav¬ 
ing  the  same  <uc  and  <up  and  driven  with  the  same  input 
currents.  For  example,  if  the  single  junction  shows 
chaotic  behavior,  so  will  the  array  and  a  coupled  state 
will  result. 

It  is  worth  commenting  on  the  behavior  should  the  in¬ 
put  currents  /,( j)  not  be  independent  of  j,  which  is  likely 
to  be  the  case  for  a  real  tray.  We  have  done  studies  with 
the  individual  currents  i ,(_/')  varied  by  random  amounts 
5/,(  j)  relative  to  a  given  mean  value  <,0,  employing  6 /',( j)  j 

as  large  as  0.5.  Then  the  currents  in  the  array  cannot  ful¬ 
ly  decouple  into  a  set  of  parallel  line  currents.  Rather, 
states  which  would  be  decoupled  for  uniform  inputs  tend 
to  show  transverse  currents  on  the  order  of  the  variation  1 
6 <,(  j )  at  the  edges  where  the  input  and  output  take  place;  ; 

these  transverse  currents  decay  toward  the  center  of  the  : 

array.  We  find  this  to  be  the  case  even  for  large  8 it(j), 
relative  to  the  distance  from  a  boundary  between  the 
decoupled  and  coupled  dynamical  states,  indicating  that 
the  qualitative  response  of  the  array  is  quite  stable 
against  nonuniformities  in  the  inputs.  An  example  is  pro¬ 
vided  by  arrays  corresponding  to  Fig.  2(c)  with  <^  =  1.0. 
using  <,o  =  0.32,  which  is  less  than  0.02  above  a  region  of 
coupled  dynamical  states,  we  find  that,  even  with  varia¬ 
tions  of  the  individual  input  currents  as  large  as  0.5,  the 
dynamics  remains  effectively  decoupled  and  periodic  with 
the  transverse  currents  becoming  small  in  the  center  of 
the  array;  for  an  8X8  array,  these  currents  are  already 
smaller  at  the  center  than  at  the  edges  by  an  order  of 
magnitude.  We  have  studied  this  effect  by  computing  the 
mean  value  of  the  absolute  current  through  transverse 
junctions  in  larger  arrays  and  find  that  it  decreases  ex¬ 
ponentially  with  the  distance  from  the  end  (<  =  1 )  of  the  j 

array.  For  an  array  with  M  =8,  the  decay  length  is  on  j 

the  order  of  2  4  junction  spacings.  This  length  is,  some-  j 

what  surprisingly,  quite  insensitive  to  /?,.  (3r.  <,0,  or  <,x.  It  I 

is,  in  fact,  a  size  effect  and  is  proportional  to  the  array  [ 
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width  M.  Hence,  our  conclusion  is  that,  for  a  sufficiently 
long  (large-AO  array,  the  global  character  of  the  long-time 
stable  state  is  quite  insensitive  to  nonuniformities  in  the 
dc  bias  current;  rather,  it  depends,  very  sensitively,  on 
the  total  input  current. 

We  give  next  a  more  detailed  description  of  the  charac¬ 
ter  of  the  different  possible  types  of  dynamical  states, 
coupled  and  decoupled,  including  representative  exam¬ 
ples  of  each  kind  of  behavior.  The  decoupled  states,  by 
current  conservation,  must  be  such  that  the  currents  in 
all  of  the  longitudinal  junctions  are  the  same  at  all  times. 
If  one  looks  separately  at  the  supercurrent  and  normal 
current  or  voltage,  however,  one  finds  that  in  some  in¬ 
stances  these  are  the  same  in  all  junctions  in  the  line  at  a 
given  time  and  in  other  instances  they  are  not.  In  the 
former  case,  the  array  is  simply  a  coherent  superposition 


of  single  junctions  and  we  refer  to  the  dynamical  state  as 
a  “single-junction”  state;  it  is  also  known16  as  a  phase- 
locked  state.  In  the  latter  case,  the  supercurrents  and 
voltages  in  different  junctions  in  the  line  are  the  same,  but 
they  are  shifted  in  time  relative  to  one  another.  We  have 
given  the  name  “chain  states”  to  dynamical  states  of  this 
kind;  they  are  also  found  in  one-dimensional  arrays 
where  they  are  known16  as  antiphase  solutions.  Note, 
however,  that  the  existence  of  such  states  in  two- 
dimensional  arrays  is,  at  first  sight,  quite  remarkable;  in 
spite  of  the  fact  that  the  phase  differences  across  two  ad¬ 
jacent  junctions  in  a  longitudinal  chain  do  not  vary  with 
time  in  the  same  way,  no  current  or  phase  difference  is 
produced  across  any  transverse  junction. 

For  the  single-junction  or  phase-locked  states,  we  find, 
for  all  parameter  sets  investigated,  that  the  period  is 


1  1  s  in  us  on  us  |  o  1  s 

sin  o' 

FIG  3  For  a  single-junction  state  at  (5,  —  4,  (},  ~  iri/25.  =  1.0,  and  i,  =  0.4,  we  show  for  a  longitudinal  junction  (a)  the  super- 

urrent  as  a  function  of  r,,  'hi  the  power  spectrum  of  ri>',  and  'cl  the  phase-space  trajectory  <  <5 '  vs  All  longitudinal  junctions  are 
I'lnalent  in  this  state. 
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equal  to  the  period  of  rf  field  or  to  some  integral  multiple  single  junction17  with  the  same  parameters,  except,  of 

of  it.  As  an  example  of  this  type  of  behavior,  we  present  course,  the  chaotic  states  involve  a  single  junction, 

in  Fig.  3(a)  the  supercurrent,  Fig.  3(b)  the  Fourier  trans-  As  an  example  of  a  chain  state,  we  present  in  Fig.  4(a) 
form  or  power  spectrum  of  4> ',  and  Fig.  3(c)  the  phase-  the  supercurrent,  Fig.  4(b)  the  power  spectrum  of  <$ ',  and 

space  trajectory  <b '  versus  sin <f>'  for  a  longitudinal  junc-  Fig.  4(c)  the  phase-space  trajectory  for  a  longitudinal 

tion  in  an  array  with  (3C  =4,  0r  =  ir /25,  ie  =0.4,  and  junction  in  the  interior  of  an  array  with  (3C  =4, 

irx  =  1.0;  cf..  Fig.  2(c).  There  is  no  mistaking  the  simple  /3f  =  ir2/25,  it  =  1.8,  and  i„  =  1.0.  Insofar  as  we  can  tell, 

periodic  character  of  the  motion  with  a  period  of  2n/(or.  the  motion  is  not  precisely  periodic.  If  it  is  periodic,  the 

If  one  increases  ie,  leaving  other  parameters  the  same,  period  is  much  larger  than  the  period  of  the  driving  field, 

period  doubling  is  observed  at  ie» 0.59,  and  again  at  This  is  a  characteristic  property  of  chain  states  in  gen- 

0.64;  there  is  further  period  doubling  followed  by  the  on-  eral.  The  power  spectrum  shows  large  apparently  sharp 

set  of  chaotic  behavior — and  a  coupled  dynamical  peaks  at  frequencies  co/<jr  equal  to  1.704,  2.704,  etc.;  the 

state — around  ie=  0.66.  Similarly,  if  one  decreases  ie,  a  widths  of  these  peaks  are  smaller  than  the  resolution  of 

coupled  dynamical  state  is  found  to  appear  by  the  time  ie  our  calculation.  There  are  smaller  sharp  peaks  at  integral 

has  reached  0.30.  Much  the  same  behavior  appears  in  a  values  of  co/a)r.  Finally,  there  are  much  smaller  peaks  at. 


~,  0|  5  -|()  -0  5  OO  0  5  II)  15 

sin  0' 

F  IG.  4  For  a  chain  state  at  (3,  =4,  /?,  -  rr/ 25,  #„  =  1.0,  and  ir  =  1.8,  we  show  for  the  longitudinal  junction  between  nodes  (4,4)  and 
1 5.4i  <ai  the  supercurrent  as  a  function  of  r,,  (b)  the  power  spectrum  of  i>',  and  (c)  the  phase-space  trajectory.  The  last  of  these,  as 
sv ell  as  all  phase-space  trajectories  in  the  following  figures,  is  discrete,  the  points  being  at  intervals  of  the  step  size  used  when  integrat¬ 
ing  the  equations  of  motion 
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e.g.,  (1—0.704),  (1+2X0.704),  and  (1  +  3X0.704),  that  the  supercurrents  or  normal  currents  in  different 

which  suggests  that  there  are  two  basic  frequencies  of  junctions  along  a  longitudinal  chain  are  the  same  but 

motion,  one  of  which  is  certainly  m,  and  the  other  of  shifted  in  time  from  one  another.  We  have  found,  by 

which  is  around  0.704cu,.  The  system  generates  frequcn-  starting  the  system  with  judiciously  chosen  initial  condi- 

cy  components  which  are  various  linear  combinations  of  tions,  that  locally  stable  dynamical  states  are  obtained  in 

these  basic  frequencies.  If  the  two  frequencies  are  incom-  which  this  time  shift  is  a  multiple  of  the  driving  period 

mensurate,  then  the  motion  is  not  truly  periodic.  Similar  2n/cjr.  We  conjecture  that  all  locally  stable  chain  states 

behavior  is  observed  for  other  chain  states;  that  is,  the  have  this  property  but  have  not  been  able  to  verify  our 

power  spectrum  has  sharp  peaks  at  integral  values  of  conjecture  in  the  case  of  dynamical  states  obtained  from 

a)/(o,  and  also  at  one  or  more  sets  of  additional  frequen-  random  initial  conditions.  In  any  event,  the  basin  of  at- 

cies  displaced  from  the  former  set  by  fixed  amounts  traction  for  any  given  chain  state  is  small  in  that  the  time 

which  do  not  appear,  in  general,  to  be  commensurate  shifts  of  the  different  junctions  depends  sensitively  on  the 

with  the  driving  frequency  &>,.  Furthermore,  as  the  junc-  initial  conditions.  The  basic  motion  of  any  given  junc¬ 
tions'  parameters  and  input  currents  are  varied,  the  tion,  however,  does  not.  Also,  an  individual  chain  state  is 

second  basic  frequency  also  varies,  apparently  smoothly.  typically  easily  disrupted  with  a  small  perturbation  to  the 

As  stated  above,  a  general  feature  of  the  chain  states  is  system  in  that  the  time  shifts  change  as  a  result  of  such  a 
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FIO.  5  For  a  quasipenodic  coupled  state  at  /?,  =  I,  0.  =  25,  /„  =0.5.  and  i,  =  1.25,  we  show  (a)  the  supercurrent  as  a  function  of  t, 
i  'he  longitudinal  junction  between  nodes  <4,4i  and  t5.4i.  'hi  the  power  spectrum  of  <j>'  in  the  same  junction,  (c)  the  phase-space  tra- 
-•orv  for  this  junction,  and  ' d i  the  phase-space  trajectory  for  the  transverse  junction  between  nodes  (4,4)  and  (4,5). 
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F  IG  6.  For  a  chaotic  state  at  0,  =4,  0,  =  n2/25,  i„  =  1.0,  and  i,  =0.69,  we  show  (a)  the  supercurrent  as  a  function  of  r,  in  the  lon¬ 
gitudinal  junction  hetween  nodes  (4,4)  and  (5,4).  lb)  the  power  spectrum  of  i>‘  in  the  same  junction,  (c)  the  phase-space  trajectory  for 
this  junction,  (di  the  power  spectrum  of  the  neighboring  transverse  junction  between  nodes  (4,4)  and  (5,4),  and  (e)  the  phase-space  tra- 
lectory  for  that  transverse  junction. 
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perturbation.  Further  properties  of  these  states,  such  as 
their  stability,  will  be  addressed  in  a  separate  publication. 

As  for  the  coupled  states,  we  have  examined  the 
motions  for  many  different  parameter  sets  and  have  seen 
just  three  qualitatively  different  types  of  long-time  behav¬ 
ior,  these  being  periodic,  quasiperiodic,  and  chaotic.  The 
first  of  these  have  periods  which  are  integral  multiples  of 
ln/a)r.  For  example,  at  /3C=4,  /3r  =  rr1/25,  /,=  1.235, 
and  i „  =  1.0,  there  is  a  coupled  periodic  state  with  period 
5(2) T/cor).  These  periodic  states  also  have,  in  general, 
some  spatial  periodicity  which  is  typically  two  junction 
spacings.  In  any  given  state  this  pattern  may  be  repeated 
throughout  the  array  or  may  be  only  local  with  different 
periodic  (in  space)  states  present  in  different  parts  of  the 
array.  Finally,  the  coupled  periodic  state  tend  not  to  be 
very  robust,  meaning  that  they  usually  have  relatively 
small  basins  of  attraction  and  are  not  found  for  most  sets 
of  initial  conditions.  Rather,  in  those  places  where  these 
states  can  appear,  there  is  usually  a  more  robust  periodic 
decoupled  state. 

An  example  of  a  coupled  quasiperiodic  state  is  provid¬ 
ed  by  an  array  with  /?r  =  1,  /?,  =25,  it  =  1.25,  and 

=0.5.  Figure  5(a)  shows  the  supercurrent  in  a  typical 
longitudinal  junction  in  the  interior  of  the  array;  Fig.  5(b) 
displays  the  power  spectrum  of  <j> '  for  the  same  junction; 
and  Figs.  5(c)  and  5(d)  present,  respectively,  the  phase- 
space  trajectories  for  this  junction  and  for  an  adjacent 
transverse  junction.  In  addition  to  peaks  at  integral  mul¬ 
tiples  of  the  driving  frequency,  there  is  a  pronounced 
peak  at  a  frequency  close  to,  but  not  precisely  equal  to, 
cjp  =w{.  =0.2oj.  Further,  there  are  smaller  peaks  at  har¬ 
monics  of  this  one's  frequency.  These  peaks  are  narrower 
than  the  resolution  of  our  Fourier  transform.  As  was  the 
case  for  the  chain  states,  the  motion  appears  to  be  quasi- 
periodic  rather  than  truly  periodic. 

A  clear  contrast  to  Fig.  5  is  shown  by  Fig.  6  where  we 
plot  in  Ftg.  6(a)  the  supercurrent  in  a  longitudinal  junc¬ 
tion,  Fig.  6(bi  the  power  spectrum  of  6'  in  that  junction. 
Fig.  6(c)  the  phase-space  trajectory  for  the  same  junction, 
Fig.  6<d)  the  power  spectrum  of  a  neighboring  transverse 
junction,  and  Fig.  6(e)  the  phase-space  trajectory  for  that 
transverse  junction,  in  an  array  with  /?r  =  4,  /J,  =  tr2/ 25, 
/. ,=0.69,  and  irx  =  1.0.  This  is  a  chaotic  state  as  is 
demonstrated  by  the  power  spectrum  and  the  phase-space 
trajectories.  The  power  spectrum  also  shows  sharp  peaks 
at  harmonics  of  the  drivirj  frequency;  these  are  absent 
from  the  power  spectrum  for  a  transverse  junction. 

For  applications,  one  is  interested  not  so  much  in  the 
behavior  of  individual  junctions  in  the  array,  but  rather 
in  global  properties,  in  particular,  the  current-voltage 
characteristics.  In  Fig.  7(a)  we  display  the  time-averaged 
mean  voltage  across  an  8X8  array  as  a  function  of  /,.. 
The  parameters  for  this  system  are  fl .  = 4,  pr  =  n2/25. 
and  i„  =  1.0.  The  curve  shows  a  number  of  voltage  pla- 
'eaus  or  Shapiro  steps  where  the  mean  voltage  is  indepen¬ 
dent  of  the  input  current.  In  these  regions  the  dynamical 
'-tales  are  the  decoupled  states  with  single-junction  be¬ 
havior;  compare  Fig.  2<c).  There  is  a  wide  plateau  for 
'*  315"  i  10.65  and  smaller  ones  in  other  regions  includ¬ 
'd.;  many  at  small  i,  ,  as  in  the  case  of  a  single  junction,11 
•vhich  are  two  narrow  to  appear  in  the  figure.  There  are 


no  plateaus  for  it  larger  than  about  1.5,  although  this  is  a 
region  where  the  dynamical  states  are  decoupled.  These 
states  are,  however,  chain  states  rather  than  single¬ 
junction  states.  Floquet  analysis,  similar  to  that  done  in 
Ref.  10,  demonstrates  that  the  single-junction  states  are 
unstable  here.  As  a  rule,  chain  states  correspond  to  re¬ 
gions  where  the  voltage  varies  with  ie.  Figure  7(b)  shows, 
for  contrast,  the  I-  V  curve  for  a  single  junction  with  the 
same  parameters  as  used  for  Fig.  7(a).  The  structure  of 
the  curve  is  similar  but  with  more  pronounced  voltage 
variations.  Also,  of  course,  the  voltage  is  seven  times 
larger  for  the  array. 

These  I-  V  curves  may  be  contrasted  with  what  is  found 
if  the  capacitance  is  zero,  In  that  case,  there  is  a 

sequence  of  wide,  evenly  spaced  voltage  plateaus  for 
currents  up  through  and  well  beyond  what  is  shown  in 
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FIG.  7.  For  (3t  =4,  (3,  -  ir/25,  and  =  1.0,  we  show  the  l-V 
curves  for  (a)  an  8X8  array  and  (b)  a  single  junction.  The  volt¬ 
ages  on  corresponding  plateaus  are  seven  times  larger  for  the  ar¬ 
ray  than  for  the  single  junction. 
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the  figures.  Hence,  we  can  conclude  that  the  addition  of 
capacitance  to  the  junctions  in  this  instance  suppresses 
the  higher  Shapiro  steps,  at  least  for  the  dynamical  be¬ 
havior  at  the  times  which  we  investigated.  At  the  same 
time,  the  edges  of  the  remaining  steps  tended  to  become 
sharper  with  the  addition  of  the  capacitance,  reflecting 
the  fact  that  the  presence  of  capacitance  does  not,  in  gen¬ 
eral,  result  in  an  overall  degradation  of  the  steps,  as  has 
been  demonstrated  also  by  others.12,18  In  the  present  ex¬ 
ample,  as  described  above,  there  is  just  one  wide  plateau 
remaining  with  a  number  of  narrower  ones  which  do  not 
necessarily  occur  at  voltages  that  are  integral  multiples  or 
simple  fractions  of  the  voltage  of  the  wide  one  (which  is 
at  the  same  voltage  as  the  first  step  produced  by  an  array 
with  negligible  capacitance). 

In  the  regions  between  plateaus,  where  the  voltage  typ¬ 
ically  increases  with  the  dynamical  state  is  either  a 
chain  state  or  a  coupled  state.  Specifically,  at  relatively 
large  ie,  ie~1.5  and  above,  it  is  the  chain  state  that  is 
present  while  for  ie  in  the  region  below  about  1.0,  there 
are  chaotic  states.  Thus  there  is  a  simple  correspondence 
of  decoupled,  periodic  dynamical  states  with  the  voltage 
plateaus  and  either  coupled  states  or  aperiodic  decoupled 
states  with  the  regions  where  the  voltage  changes  with  in¬ 
put  current.  Similar  statements  apply  in  this  case  of  a 
single  junction;  there  are  periodic  states  on  the  voltage 
plateaus  and  either  chaotic  or  quasiperiodic  states  else¬ 
where. 

We  have  studied  the  transition  from  periodic  to  chaot¬ 
ic  behavior  in  some  detail.  As  an  example,  we  cite  the 
transitions  as  one  moves  from  the  large  plateau  in  Fig. 
7(a)  either  upward  or  downward  in  ie  into  chaotic  re¬ 
gimes.  In  the  former  case,  there  is  a  sequence  of  period 
doubling  transitions  with  increasing  ie  until  chaotic  be¬ 
havior  sets  in  around  /,  =0.66,  a  standard  route  to  chaos. 
The  behavior  for  decreasing  /t.  is  not  the  same.  In  this 
case  we  find  no  period  doubling.  Rather,  the  power  spec¬ 
trum,  w  hich  is  necessarily  initially  a  set  of  perfectly  sharp 
peaks  at  integral  multiples  of  <or,  evolves  by  developing 
widths  in  the  peaks  accompanied  by  the  gradual  appear¬ 
ance  of  a  noisy  background  at  all  frequencies.  Hence, 
two  quite  different  routes  to  chaos  appear  in  this  one  ex¬ 
ample. 

It  is  appropriate  to  compare  the  behavior  of  the  array 
with  that  of  a  single  junction  having  the  same  charac¬ 
teristics.  In  many  respects  the  two  are  comparable.  That 
is.  the  single  junction  also  displays  periodic,  quasiperiod¬ 
ic,  and  chaotic  behavior.  Of  course,  with  a  single  junc¬ 
tion  one  cannot  observe  coupled  as  opposed  to  decoupled 
dynamical  states.  Roughly  speaking,  the  regimes  where 
the  array  shows  coupled  chaotic  states  are  also  ones 
where  the  single  junction  is  chaotic;  where  the  array 
displays  periodic  behavior,  the  single  junction  is  also 
periodic  with  the  same  period,  and  where  the  array  shows 
aperiodic  but  decoupled  behavior,  the  single  junction  is 
also  aperiodic.  It  is.  however,  sometimes  difficult  to 
make  a  detailed  comparison.  As  a  rule,  if  the  single  junc¬ 
tion  is  chaotic  then  the  array  is;  there  are  some  regions 
where  the  array  is  apparently  chaotic  but  the  single  junc- 
1 1<  >ii  is  not  This  observation  may  be  a  consequence  of  the 
finite  limes  we  are  able  to  follow  the  motions  of  the  array. 


That  is,  it  is  not  inconceivable  that  in  some  of  these  cases 
the  array  would  exhibit  long-time  periodic  or  quasi¬ 
periodic  behavior  were  we  able  to  follow  it  long  enough. 
Consider  as  an  example  junctions  with  0C  =  4,  0r  =  it1/25, 
and  i„  =  1.0.  For  ie  between  0.07  and  0.13,  there  is  a  re¬ 
gion,  corresponding  to  the  voltage  plateau  in  Fig.  7(b), 
where  the  single  junction  is  periodic.  At  many  values  of 
ie  in  this  range  we  found,  for  8  X  8  arrays,  chaotic  rather 
than  periodic  behavior  starting  from  randomly  selected 
initial  states. 

B.  Effects  of  noise 

The  addition  of  noise  currents  to  the  system  provides  a 
means  by  which  it  may  escape  from  one  locally  stable 
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sin  0' 

FIG.  8.  The  phase-space  trajectory  for  (a)  the  longitudinal 
junction  between  nodes  (4,4)  and  (5,4)  and  (b)  the  transverse 
current  in  the  junction  between  nodes  (4.4)  and  (4.5),  for  =4. 
(3.  =  -:/25.  i.t  -  1.0.  i  =0.32.  and  1 1, .  50.4. 
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dynamical  state  to  another.  Further,  the  characters  of 
various  types  of  states  must  change  somewhat; 
specifically,  the  decoupled  states  are  no  longer  truly 
decoupled  because  fluctuating  currents  must  appear  in 
the  transverse  junctions.  In  addition,  periodic  states  are 
no  longer  truly  periodic  because  of  the  noise.  The  most 
interesting  question  associated  with  the  presence  of  noise 
is  whether  states  are  destabilized  by  it. 

One  well-known19  consequence  of  the  noise  is  that  the 
edges  of  the  plater  *:s  in  the  current-voltage  plot  become 
rounded.  Small  plateaus  are  lost  altogether.  Purely 
periodic  states  become  quasiperiodic  in  that  the  power 
spectrum  still  has  sharp  peaks  with  some  background 
present.  The  phase-space  trajectories  demonstrate  the 
quasiperiodic  behavior  very  vividly;  Fig.  8,  which  may  be 
compared  with  Fig.  3(c),  is  for  /3C  =4,  (3r  =  1^/25, 
4,  =  1.0,  and  4=0.32,  with  the  rather  large  noise 
current:o  |/'f(f)|  50.4.  Even  this  large  fluctuating  current 
is  not  sufficient  to  disrupt  the  basic  periodic  state;  notice, 
however,  that  in  the  absence  of  noise,  if  ie  were  reduced 
to  0.30,  a  reduction  of  just  0.02,  a  chaotic  state  would  re¬ 
sult.  The  periodic  unperturbed  state  is  evidently  quite 
robust  against  noise,  even  as  it  is  robust  against  nonuni¬ 
form  input  currents. 

Comparably  large  noise  currents  also  do  not  destroy 
other  types  of  unperturbed  states.  We  mention,  in  partic¬ 


ular,  the  chain  states  which  maintain  their  basic  charac¬ 
ter,  described  above,  in  the  presence  of  the  noise. 


IV.  STABILITY  ANALYSIS 

From  Sec.  II,  the  condition  of  current  conservation  at 
node  k  may  be  expressed  as,  for  iF= 0, 

2 

/6vU) 

+sin(^-#)+i;(/c)=0  .  (27) 


i/2  d 
dr 


In  the  following,  we  shall  specify  a  node  using  indices  ij 
as  in  Sec.  II,  <p'k(r)— ;r). 

We  begin  by  considering  the  stability  of  the  single¬ 
junction  states.  These  are  characterized  by  a  single-phase 
function  <£0(r).  If  we  let  the  phases  of  the  individual 
nodes  in  the  single-junction  state  be  designated  as 
<j>'0{ij;r),  then  these  are  related  by 


<P'Q( i,j ; r) — <t>'Q( i',j';T)  =  ( /  —  / ' )<p0(r)  .  (28) 

To  test  this  mode’s  stability,  we  write  <p'(i,j\r)  as 

<f>'( i,j ; r)  =  <£q( ij ;T)  +  \p(i,j  ;t)  ,  (29) 


and  expand  Eq.  (27)  in  powers  of  ip,  keeping  just  linear 
terms. 


[  4t /<(  i,j  )  —  \p(i  +  l,j)  —  <pU  —  l,j)  —  \p(  i,j  +  1 )  —  ip{  i,j  —  1 )  ] 


+  cos<60(  2tM  /,  j )  - «/-( /  +  1 ,  j )  -  4>(  i  ~  I ,  j )  ]  +  [  2 ip(  i , ; )  -  / ,j  +  1 )  —  tpi  i,  j  -  1 )  ] = 0  .  (30) 


r 


Finally,  we  Fourier  transform  the  linearized  equations  of 
motion  in  space,  making  use  of  the  fact  that  the  nodes  are 
on  a  square  lattice  with  spacing  a,  and  obtain 


dr 


:  d  i/'l  q;  r ) 


-/[ q; <*r,(  r)]t/'(q;r)=0 


(31) 


*2  = 


dtpir) 

dr 


=  i'  , 


so  that 

*',=x2 


where 


and 


/ 


i  I  —cos qxa  )cos (6,,-H  1  —cos qva  ) 
2  —  cost?,  a  —  cos  qva 


(32) 


*2 - K 


1/2, 


'I*  I 


In  the  single-junction  mode,  cos[<Jn(r)]  is  periodic,  and 
as  a  consequence  /(q,r)  is  also  periodic  so  that  Floquet 
theory:l  may  be  applied.  First,  note  that,  for  1  qx ,  <  I qy  l, 
/  q, r)>0  for  all  r.  Consequently,  the  uniform  solution 
is  stable  against  perturbations  which  satisfy  this  condi¬ 
tion  on  q.  More  generally,  we  convert  the  second-order 
equation  ' 3 1 1  to  a  pair  of  first-order  equations  by  intro¬ 
ducing 

x  .  i  r  i  =  tl’1  r  i 


and 


1 33) 


or,  in  matrix  form, 


x'  —  Ax  , 


where 


x  = 


and 


0  I 

-/<->  -Pr 


(34) 


(35) 


(36) 


A  = 


1.2 
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We  have  analyzed  these  equations,  finding,  in  particu¬ 
lar,  the  two  characteristic  multipliers  a,  for  a  variety  of 
single-junction  states.  In  all  cases  examined  we  find  that 
the  absolute  values  I  a,  |  are  less  than  unity  for  any  finite 
qv.  However,  one  of  them  is  very  close  to  unity  for  qy 
small.  The  implication21  is  that  the  periodic  uniform 
states  are  (locally)  asymptotically  stable,  although  the  de¬ 
cay  time  of  a  perturbation  can  be  very  long.  These  pre¬ 
dictions  are  supported  by  our  numerical  results;  all 
periodic  uniform  states  that  we  have  encountered  have 
proved  to  be  locally  stable. 

Next,  suppose  that  cos[0o(  r )  ]  is  chaotic  corresponding 
to  a  decoupled  uniform  state  in  which  each  junction  un¬ 
dergoes  the  same  chaotic  motion.  In  this  case  the  uni¬ 
form  solution  is  one  in  which  the  entire  system  behaves 
like  a  single  junction  and  exhibits  uniform  global  chaos. 
We  may  argue,  however,  that  this  state  is  unstable  with 
respect  to  small  perturbations.  Suppose  a  phase  at  one 
node  is  perturbed  relative  to  its  value  in  the  uniform 
state.  This  disturbance  will  grow  as  a  consequence  of  the 
state  being  a  chaotic  one  with  a  positive  Liapunov  ex¬ 
ponent.  Further,  the  disturbance  will  spread  to  the  near¬ 
by  nodes,  which  are  coupled  to  the  original  one,  and  so 
will  grow  in  spatial  extent  as  time  progresses  as  well  as 
growing  on  each  node  where  it  is  present;  hence,  the  uni¬ 
form  state  is  unstable.  This  conclusion  is  also  supported 
by  our  numerical  results.  All  chaotic  states  that  we  have 
studied  show  short-ranged  spatial  correlations  in  the 
chaotic  part  of  the  motion.  However,  the  chaotic  fluctua¬ 
tions  tend  to  be  accompanied  by  some  underlying  period¬ 
ic  motion  with  longer-ranged  correlations. 

V.  DISCUSSION  AND  SUMMARY 

We  have  studied  the  dynamical  states  of  two- 
dimensional  arrays  of  capacitive  Josephson  junctions 
with  dc  bias  and  an  applied  monochromatic  rf  field  using 
the  RSJ  model.  For  wide  ranges  of  junction  parameters, 
we  find  stable  or  metastable  solutions  of  several  basic 
types.  First,  the  states  are  either  coupled,  meaning  that 
there  are  currents  in  junctions  perpendicular  to  the  input 
currents,  or  decoupled,  meaning  that  the  input  currents 
pass  straight  through  the  array,  producing  no  currents  in 
the  transverse  direction.  In  the  latter  case,  the  arrays 
display  either  single-junction  (phase-locked)  states  or 
chain  states  in  which  the  junctions  in  a  longitudinal  chain 
all  behave  in  the  same  fashion  but  successive  junctions  in 
the  chain  are  time  shifted  by  apparently  random  multi¬ 
ples  of  the  period  of  the  rf  field.  At  the  same  time,  the 
states  of  parallel  junctions  in  different  chains  are  identical 
so  that  the  system  is  invariant  under  translation  in  the 
transverse  direction. 

Also,  both  kinds  of  decoupled  states  are  generally  very 
robust  in  that  random  changes  in  the  initial  conditions  do 
not  often  lead  to  a  different  dynamical  state  at  long  times; 
in  the  case  of  the  chain  state,  such  a  perturbation  typical¬ 
ly  produces  a  different  chain  state  in  that  the  amounts  by 
which  different  junctions  are  time  shifted  relative  to  each 
other  change.  Similarly,  thermal  noise,  included  by 
means  of  a  Langevin  noise  current  in  each  junction,  does 
not  easily  destabilize  these  states  although  it  does,  of 


course,  cause  them  to  fluctuate  around  some  mean  prop¬ 
erties.  Similarly,  random  variations  in  the  input  currents 
do  not  have  significant  effects  on  the  decoupled  states  so 
long  as  the  total  input  current  is  unchanged.  In  the  latter 
case,  there  are  produced  currents  in  transverse  junctions 
close  to  the  nodes  where  the  current  is  injected,  but  these 
die  out  exponentially  with  distance  with  a  decay  length 
on  the  order  of  one-half  of  the  array’s  transverse  length. 

The  rather  impressive  stability  of  the  dynamical  modes 
against  thermal  noise  and  disorder  in  the  inputs  is  some¬ 
what  surprising.  It  would  be  useful  to  know  also  whether 
the  dynamical  states  are  as  stable  when  there  is 
configurational  disorder  in  the  array;  we  intend  to  study 
this  question  and  report  on  it  separately.  We  do  not  at 
this  point  have  an  explanation  for  this  stability.  Never¬ 
theless,  it  is  encouraging  in  that  it  implies  that  reliable 
devices  such  as  rf  radiation  detectors  could  be  construct¬ 
ed  from  junction  arrays. 

The  coupled  dynamical  states  are  either  periodic, 
quasiperiodic,  or  chaotic.  The  periodic  states  that  we 
found  to  be  uncommon  in  that  they  have  relatively  small 
basins  of  attraction  and  do  not  appear  very  often  if  the 
simulation  is  started  with  random  initial  conditions. 
More  commonly,  for  those  parameters  where  these  states 
can  appear,  a  periodic  decoupled  state  usually  is  the  re¬ 
sult  of  starting  from  an  arbitrarily  chosen  configuration. 
The  periodic  states  are  interesting  in  that  they  also 
display  spatial  periodicity.  As  for  the  quasiperiodic 
states,  these  behave  in  time  rather  like  the  chain  states 
but  have,  in  addition,  currents  in  the  transverse  junctions. 
Finally,  chaotic  states  are  quite  common;  we  have  found 
that  chaos  can  appear  by  a  sequence  of  bifurcations  dur¬ 
ing  which  the  array  is  in  a  single-junction  state  and  also 
by  intermittency. 

The  current-voltage  characteristics  of  the  array  have 
been  examined  and  correlated  with  the  dynamical  states. 
There  are  voltage  plateaus  where  the  voltage  remains 
constant  for  some  range  of  ie  or  i rx.  The  dynamical  states 
on  these  plateaus  are  single-junction  states  or  periodic 
coupled  states.  Off  the  plateaus  one  finds  coupled  and 
chain  states.  The  presence  of  capacitance  in  the  arrays 
tends  to  suppress  the  plateaus  and  produce  fairly  smooth 
variation  of  the  voltage  with  the  input  currents,  especial¬ 
ly  in  the  region  of  relatively  large  ie  (larger  than  about 
1.5)  where  predominantly  chain  states  are  found.  For 
smaller  :t.  and  with  irx  ~  1,  chaotic  states  are  not  common 
off  the  voltage  plateaus.  Also,  in  this  region  there  are 
numerous  very  small  plateaus.  Overall,  the  behavior  of 
the  I-V  curves  is  not  unlike  that  for  a  single  junction  with 
the  same  properties  as  the  junctions  in  the  array,  al¬ 
though  much  of  the  detail  of  the  curve  for  the  single 
junction  is  lost. 

Of  the  various  phenomena  we  have  found,  several  seem 
worthy  of  further  investigation.  These  include  the  de¬ 
tailed  character  of  the  chaos  in  this  extended, 
continuous-time  system,  the  rather  remarkable  stability 
of  the  various  dynamical  states,  especially  the  decoupled 
ones,  to  different  kinds  of  noise  and  disorder,  and  certain 
features  of  the  chain  states.  Further  work  on  these  and 
other  aspects  of  the  dynamics  of  capacitive  arrays  is  in 
progress. 
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